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The high-frequency peaked BL Lac object 1ES0229H-200 (z = 0.14) was first detected in very 
high energy (VHE, E>100GeV) y-rays by the H.E.S.S. (High Energy Stereoscopic System) 
collaboration in 2006 [1]. No flux variability was reported in the initial study and its spectral 
characteristics have been used to derive constraints on the extragalactic background light (EBL, 

[1]) and on the intergalactic magnetic field (IGME, e.g. [2, 3, 4, 5]). 1ES0229H-200 has been 
observed with H.E.S.S. for ~130 hours from 2004 to 2013: the full dataset analysed with a more 
sensitive method will be presented here. The results indicate that the source is not constant and 
displays flux variability on yearly and monthly timescales. The existence of flux variability affects 
the derivation of the constraints on the IGME. The H.E.S.S. observations cover several simulta¬ 
neous multi-frequency campaigns and the VHE variations are compared with those reported in 
different bands. 
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1. Introduction 

lES 0229+200 is a blazar, and a member of the high frequency peaked BLLac objects (HBL) 
subclass. It is hosted in an elliptical galaxy located at a redshift z = 0.1396. It was first discovered 
in X-rays by the Einstein IPC Slew Survey in 1992 [6]. Its synchrotron emission reaches 100 keV 
and possibly up to 200keV^ Despite being relatively bright in optical and X-rays, it is faint in 
radio and y-rays. It has been detected at VHE above 500 GeV for the first time by H.E.S.S. in 2006 
[1], and only in 2011 at high energies (HE, 100MeV<E< 100GeV) between 1 and 300GeV, 
after more than three years of Fermi observations [5]. Contrary to what one could expect from its 
blazar nature, lES 0229+200 shows little variability at all energies, even if it is known to have flux 
variafions in fhe X-ray band of a factor of ~ 2. (e.g. [7]). Hints of variability at TeV energies were 
reported on yearly and monthly timescales by the VERITAS collaboration in 2014 [8], while none 
was detected in the source in the initial H.E.S.S. study. 

lES 0229+200 has been one of the key sources for deriving constraints on the extragalactic 
background light (EBE, [1]) and on the intergalactic magnetic field (IGME, e.g. [2, 3, 4, 5]) thanks 
to the combination of its hard spectrum (Epl ~2.5, [1, 8]) reaching ~ 10 TeV and its consider¬ 
able distance (for a TeV source). The spectrum allows one to probe the near- and mid-infrared 
wavelengths of the EBE (~ 2 - 20 jum) at high optical depths (r ~ 1 - 6 for y-ray energies between 
~ 1 -10 TeV). This is unique, since only a few other sources have been observed at these energies. 
Eor the IGME studies, a multi-TeV hard spectrum in combination with high optical depths implies 
that a considerable part of the emitted TeV photons will be absorbed via pair production through 
their interaction with the EBE and re-emitted at (typically) GeV energies, causing a surplus of HE 
emission. Eower limits on the strength of the IGME can be calculated measuring the GeV emission 
and making assumptions on the intrinsic spectrum. Since reprocessing involves emission from an 
ensemble of charged particles that are spread over an extended volume by the action of the magnetic 
held, time delays are of the order of the light travel time of the excess distance that the re-emitted 
photons must cover with respect to the line of sight to the source. Time delays are also energy 
dependent, since a given magnetic held will affect more strongly the trajectories of low-energy 
particles. This holds true as long as the Compton cooling time is of the order of the gyroradius. If 
it is too fast, re-emission is almost immediate and time delay negligible. If it is too slow, this and 
not the magnetic held determines the time delay. The observed reprocessed emission corresponds 
therefore to the time-averaged original VHE emission. Constrains on the IGME hence rely on the 
assumption of a steady VHE hux on time scales of these delays. Is the VHE hux variable, average 
values can be used as long as the variability is shorter than the reprocessing time delay. 

2. Observations and Analysis 
2.1 H.E.S.S. 

H.E.S.S. observations have been carried out almost yearly between 2004 and 2013, for a total 
of 354 observation runs^. 311 of them pass the standard quality cuts for a total livetime of ~ 133 h 

'http://bat.ifc.inaf.it/100m_bat_catalog/lO0m_bat_catalog_v0.0.htm 

^The term "run" refers to a single observation, with a typical exposure of 28 minutes. 
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Table 1: Statistics of the yearly and total H.E.S.S. observations. 


Period 

Dates 

(MJD) 

Live Time 

(h) 

ON 

counts 

OFF 

counts 

excess 

Significance 

(o-) 

Flux (>580GeV) 
(10“'^ cm“^s“') 

2004-2013 

53259-56606 

132.5 

1839 

14138 

699 

18.1 

6.2+ 0.5 

2004 

53259-53317 

3.8 

47 

349 

18 

2.9 

3.5 ±1.9 

2005 

53613-53649 

16.4 

169 

1518 

42 

3.4 

2.5+ 1.0 

2006 

53967-54088 

47.1 

750 

4577 

368 

15.8 

9.0 ±0.8 

2007 

54322-54336 

10.5 

137 

1080 

47 

4.4 

5.6+1.7 

2008 

54681-54789 

12.2 

149 

1420 

31 

2.6 

4.0+1.6 

2009 

55063-55151 

19.0 

215 

1825 

63 

4.6 

5.2+1.4 

2011 

55801-55909 

7.4 

112 

918 

35 

3.6 

1.9 ±2.0 

2013 

56514-56606 

17.8 

241 

2189 

95 

6.9 

7.4+1.3 


Note: the normalization factor between the ON and OFF area is 1/12 for all years with the exception of 2013 1/15. 


(122 h when corrected for acceptance). The mean zenith angle is 45.2° and the mean offset from 
the pointing position 0.51°. The exposure was not homogeneous during these observation years, 
ranging from 3.8 h in 2004 to 47.1 h in 2006. Refined quality selection criteria allowed for a 
larger dataset for the years 2004-2006 than the one published in [9] and [1]. Part of the 2009 and 
2013 data belong to multiwavelength (MWL) campaigns. The first one was organized with XMM- 
Newton and ATOM on August 21 and 23, 2009, the second one with NuSTAR, Swift-XRI, MAGIC 
and VERITAS on October 1, 5 and 11, 2013. Data reduction has been performed using the Model 
analysis^ [10] with Standard cuts. Eor the background determination for the spectral and temporal 
analyses, the Reflected Region Background method [11] was used. The source is detected with a 
signihcance of 18.1 standard deviations (cr) for the whole dataset. The results for the total dataset 
and for every year of observation can be found in Tab. 1, as well as in Tab.2 for the 2006 and 2009 
observing periods. 

Spectral analyses have been carried out on different data subsets using the forward-folding 
technique [12] and a simple power-law (PE) spectral shape. No spectral variability is detected on 
a yearly nor monthly timescale. The lightcurves are derived adjusting the normalization of a PE 
with index E = 2.9 to the y-ray excess in every time bin. Only energies above a common threshold 
of 580GeV are considered. The monthly and yearly lightcurves are shown in the top two panels 
of Eig. 1 and discussed in Sec. 3. All results have been cross checked with a different software and 
calibration chain. 

2.2 High Energies 

HE data have been collected by the Earge Area Telescope (EAT) onboard the Fermi satel¬ 
lite starting August 4, 2008. Data up to April 27, 2015 have been analysed with the Science- 
Tool software package version v9r33pQ. Only events belonging to the ’Source’ class within 
15° from the position of lES 0229-1-200 were selected. Moreover, cuts on zenith angle (100°), 

^software version paris-0-8-24, DSTs Prod26 
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Table 2: Statistics of the 2006 and 2009 observation periods. 


Period 

Dates 

(MJD) 

Live Time 

(h) 

ON 

counts 

OFF 

counts 

excess 

Significance 

(O') 

Flux (>580 GeV) 
(10“'^ cm“^s“') 


2006 

August 

53967-53977 

13.5 

226 

1581 

94 

7.1 

8.1+ 1.6 

September 

53994-54005 

16.0 

329 

1813 

178 

11.9 

12.2+1.6 

November 

54048-54063 

12.5 

146 

813 

78 

7.8 

6.9+ 1.4 

December 

54077-54088 

5.1 

49 

370 

18 

2.9 

6.8+ 2.0 


2009 

August 

55063-55074 

9.3 

72 

815 

4 

0.5 

-0.5+ 1.5 

October 

55115-55121 

5.1 

73 

561 

26 

3.4 

7.0+ 2.9 

November 

55145-55151 

4.6 

70 

449 

33 

4.5 

16.6 + 3.8 

Note: the non 

nahzation factor p between the ON and OFF 

area is 1/12 for all periods. The observing periods are defined 


by the lunar cycle and do not always strictly match the given month, which is given for an easier identification. 


rocking angle (52°) and distance from the Sun (5°) were applied. For the binned maximum- 
likelihood spectral analysis, the instrument response functions P7REP_SOURCE_V15 were used, 
together with the publicly available standard Isotropic and Galactic diffuse emission background 
models iso_source_v05.txt and gll_iem_v05 revl. fit. 

lES 0229-1-200 is detected with a significance of 10.5 tr in the energy range 100 MeV - 500 GeV, 
compared to the 6.7 cr between 1 and 300 GeV in [5]. The low flux in the Fenni-LKT band does 
not permit detailed temporal studies on monthly timescales. Nonetheless, the variability index^ for 
the monthly lightcurve has a value of 157 for 81 degrees of freedom, which indicates variability at 
the 5 cr level. 

2.3 X-rays 

In this work. X-ray data from Swift, RXTE and XMM-Newton are used^ (Fig. 1, third and 
fourth panels). The 30-days bin Swift-BAI lightcurve between 15 and 85keV is derived from the 
15-days bin one in the 66 months Palermo BAT Catalog^. The 2008 and 2009 XMM-Newton and 
S'wi/t-XRT, as well as the 2010 RXTE data are taken from [7], while the whole 2010-2011 RXTE 
lightcurve is taken from [13]. Finally, the most recent Swift-XRI data have been collected between 
the nights of October 1-2 and 10-11, 2013 during the aforementioned MWL campaign. 

All the available Swift-XBH data collected between 2008 and 2015 (ObsIDs 00031249001- 
00031249050 and 00080245001-00080245006) were analysed using the HEASoft software pack¬ 
age V. 6.16^ with CAEDB v. 20140120. All the events were cleaned and calibrated using the 

^TSvar = 2Z,[logT;(G') - logXi(Fconst)] where Fconst and F/ are the average source flux and the flux of the t-th time 
bin. No limits on the bin significance are required. 

^NuSTAR data are not taken into account here since they are subject of a dedicated paper, in preparation. 

®http: / /bat. ifc. inaf. i t/bat_catalog_web/66iti_bat_catalog. html 
^http://heasarc.gsfc.nasa.gov/docs/software/lheasoft 


4 















Long term TeV lightcurve of lES 0229+200 


Gabriele Cologna 


xrtpipeline task and the data in the 0.3-10keV energy range with grades 0-2 for WT mode 
and 0-12 for PC mode were analysed. The lightcurve flux points were calculated from the spectra 
of single snapshots integrating between 2 and lOkeV. These were derived in the following way: 
the data were grouped using the grappha tool to have a minimum of 30 counts/bin and then fit 
using XSPEC v. 12.8.2 with a single power-law model and Galactic hydrogen absorption fixed to 
?i//-8.06xl02°cm“2 [14]. 

2.4 Optical 

Optical monitoring of lES 0229-1-200 has been carried out with the ATOM telescope [15] 
between 2007 and 2012. ATOM is a 75 cm altazimuth telescope equipped with an Apogee Alta 
U47 camera (Apogee Alta E47-I- between 2004 and 2011) and works in robotic mode. It is part 
of the H.E.S.S. project and is used mainly for optical monitoring of variable y-ray sources as a 
constant support for MWE observations. The observations for lES 0229-1-200 have been carried 
out typically once every three nights during the August-December visibility window. The standard 
reduction, the photometry and the source calibration are done with an automatic pipeline. The 
photometry uses an aperture of 4 ", while two reference stars in the vicinity of the source and present 
on the same frame are used for the calibration. During the 2013 MWE campaign ATOM was 
not operational and data were collected with the 70 cm telescope of the Abastunami Observatory 
(Georgia), equipped with an Apogee 6E camera. Observations taken with a R Cousins filter have 
been analysed with the Daophot 11 reduction software using an aperture diameter of 10 . 

3. Results and discussion 

The Model analysis is more sophisticated and sensitive than the Hillas analysis [16] used in the 
previous publications. Its use permits therefore to have a much higher sensitivity to small flux vari¬ 
ations. The monthly and yearly VHE lightcurves (top two panels of Eig. 1) show an emission which 
is not constant, neither on a yearly, nor on a monthly timescale. A fit to a constant value yields a 
of 84.3 (33.3) for 22 (7) degrees of freedom for the monthly (yearly) lightcurve. This translates 
in a probability of 5xl0~^ and 2x10“^ respectively. The values for the fractional variability ([17], 
though whose errors are calculated as in [18]) are then 0.60±0.15 and 0.39±0.11. Variability 
is detected on a monthly timescale also between 0.58 to 1 TeV and above 1 TeV (probability of a 
constant flux of lxl0~^ and 2xl0~^, respectively). 

In 2009, H.E.S.S. and VERITAS observed the source contemporaneously (though on different 
nights). The time range spanned by the VERITAS observations was significantly longer than the 
H.E.S.S. one. Taking the windows into account, the measurements of both experiments are consis¬ 
tent with an increasing flux that reaches its maximum just before MJD 55150 and then decreases. 
The comparison of the VHE and X-ray data in 2009 [7] is very interesting. The source is in very 
low state for both H.E.S.S. and XMM-Newton during the MWE campaign (MJD ~55065). This 
changes in the following three months, when Swift-XRT detects a factor ~2 flux increase [7] (also 
seen by Swift-BPLT), which is mirrored by an enhancement in the emission in y-rays of a factor 2 
or 3 (Eig. 1). In all three wavebands, this state is followed by a decrease in the emission. Despite 
the limited sensitivity of Swift-BKT in the hard X-ray band, it is interesting to notice that the three 
periods of low VHE fluxes in 2008 and 2009 are min'ored also in this band. The R and B bands 
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Figure 1: Lightcurves of the BLLac object lES 0229+200 in different energy bands between 2004 
and 2013. From top to bottom-, monthly and yearly H.E.S.S. lightcurves above 580GeV - the 
VERITAS values from [8] are also depicted as comparison; hard X-rays monthly lightcurve be¬ 
tween 15 and 85 keV from Swift-BPLT (the Palermo BAT Catalogue); soft X-ray lightcurve between 
2 and lOkeV for different instruments: 2010-2012 RXTE dataset from [13], XMM-Newton, 
Swift-XRT and RXTE data from [7], Sw//t-XRT data analysed in this work. The source is clearly 
variable in all energy bands. The discrepancy between the H.E.S.S. and VERITAS points in 2009 is 
explained by the non identical observation windows of the two instruments (see text). Three point 
at MJD 53831, 55331 and 55377 have been removed for clarity from the Swift-BKT plot because 
of their large negative fluxes or error bars. 


opfical lighfcurves, insfead, are consfanf wifhin fhe errors and do nol show any variafions. During 
fhe second MWE campaign in 2013 (MJD ~56570), lES 0229+200 seems fo be relatively brighf 
in X-rays and on ifs average value in TeV. 

A correlation in fhe X-ray-VHE emission is acfually expected by fhe SSC model, implying fhaf 
fhe fwo componenfs are generafed by fhe same elecfron populafion. On fhe ofher hand, fhe big peak 
shown by RXTE in 2011 (MJD ~55800) does nol have a clear correspondence in TeV. This could 
be explained by fhe very low exposure in Ihree oul of fhe four periods of H.E.S.S. observafions, or 
by a period of high aclivily of one of fhe ofher X-ray sources presenf in fhe field of view of RXTE 
(e.g. IRXS J023558.0+201215, IRXS J023427.5+192247 or one of fhe numerous XMM sources). 
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A third possibility could be the presence of a second X-ray emitting zone, not related to any y-ray 
emission. Finally, also Klein-Nishina suppression could play a role under certain conditions. 

The confirmation of VHE variability on these timescales (already suggested in [8]) affects the 
determination of constrains on the IGMF. Focusing on the detected variability, if the ~600 GeV 
emission consists of secondary photons ("maximal" case in [4]), then a measurement of the IGMF 
can actually be done. Extrapolating Eig. 2 of [4] to shorter time delays, one measures the IGME to 
be Biqmf ~ 3 X 10~'^ G. 

Although it has been stated that TeV y-rays from lES 0229-1-200 could be explained by proton 
cascade emission (e.g. [19]), the variability above 1 TeV implies that the radiation of this source 
is dominated by primary y-rays. Euture studies investigating further the energy dependendent 
variability will allow the contribution of the cascade emission to the flux to be probed more deeply 
for bofh hadronic and lepfonic models. 

4. Conclusion 

The 10 years long H.E.S.S. monitoring of fhe blazar lES 0229-1-200 befween 2004 and 2013 
has been presenfed for fhe firsf fime in a MWE confexf. Clear variabilify is defected af VHE on 
monfhly and yearly timescales. A hinf of correlafion befween TeV and X-ray emission comes from 
fhe confemporaneous observafions wifh XMM-Newton and wifh Swift in 2009: fhe fluxes increase 
in all energy bands in similar way over four monfhs. This supporfs an SSC emission model for fhe 
X-ray and VHE emission in fhis source. The VHE monfhly flux variabilify affecfs fhe derivation 
of lower limits on the IGME. Assuming the emission at ~600 GeV to be of secondary photons, one 
can actually obtain a measurement of the IGME of Bigmf ~ 3 x 10“'^ G. 
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